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Abstract 



The deposition of a metal on a foreign substrate is studied by means 
of grand canonical Monte Carlo simulations and a lattice-gas model 
with pair potential interactions between nearest neighbors. The influ- 
ence of temperature and surface defects on adsorption isotherms and 
differential heat of adsorption is considered. The general trends can 
be explained in terms of the relative interactions between adsorbate 
atoms and substrate atoms. The systems Ag/Au(100), Ag/Pt(100), 
Au/Ag(100) and Pt/Ag(W0) are analyzed as examples. 

Keywords: metal deposition, lattice-gas model, Monte Carlo simu- 
lation. 

1 Introduction 

The study of the adsorption of particles on surfaces is important from the 
point of view of surface science and due to its potential applications in nan- 
otechnology and catalysis. 

From the electrochemical point of view, it is of great interest to study the 
electrodeposition of a metal onto a single crystalline surface of a foreign metal 
[TJ El El H] . When this occurs at potentials more positive than those predicted 
from the Nernst equation, the process is denominated underpotential depo- 
sition (UPD) pQ El El El [7] . This can be intuitively understood considering 
that, in general, UPD is observed when for the adsorbate it is more favorable 
to be deposited on the considered substrate than on a surface of its same 
nature. When the opposite occurs, the observation of another phenomenon 
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called overpotential deposition (OPD) is expected , which implies that the 
adsorption of the adsorbate on the foreign substrate is less favorable than 
that predicted by the Nernst equation. 

It is clear that a complete analysis of the UPD phenomenon is a quite 
difficult subject because of the complexity of the involved systems. For this 
reason, the understanding of simple models with increasing complexity might 
be a help and a guide to establish a general framework for the study of 
this kind of systems. Pioneering work using Monte Carlo simulations to 
study underpotential multisite adsorption has been undertaken by Van Der 
Eerden et al. [8]. In this context, the present work tries to contribute to 
the comprehension of some essential characteristics of the metal deposition 
by means of a very simple model. For this purpose, we present simulations 
of metal adsorption on metal surfaces, employing a lattice-gas model with 
pair potential interactions characterized by a few parameters. Depending on 
the particular values assigned to the parameters, we can represent different 
metals involved. As illustrative examples we have simulated the following 
systems: Ag / 'Au(lOO) (that means: adsorption of silver atoms on gold (100) 
surfaces), Ag/Pt(100), Au/Ag(100) and Pt/Ag(100). 

In previous works P, [10] we have used a more complicated model for 
the calculation of energies, employing the embedded atom method, which 
is more appropriate for metals in the sense that takes into account many- 
body interactions. The cost of introducing this precursor model is the lack 
of some experimental features presented by real systems. However, there are 
important physical motivations to pay this cost: a) this contribution allows 
to identify and characterize the most prominent features of this particular 
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process; b) the study offers a general theoretical framework in which this 
kind of process can be studied and c) the conclusions obtained are interesting. 
The outline of the paper is as follows: In Section 2 we describe the lattice-gas 
model along with the simulation scheme. In Section 3 we present the results. 
Finally, the general conclusions are given in Section 4. 

2 Model and simulation method 
2.1 Lattice model 

Lattice models for computer simulations are of widespread use in studies of 
adsorption on surfaces, because they allow dealing with a large number of 
particles at a relatively low computational cost [TT] . 

If the crystallographic misfit between the involved atoms is not important, 
it is a good approximation to assume that the adatoms adsorb on defined 
discrete sites on the surface, given by the positions of the substrate atoms. 
This is the case of the very well studied system Ag/Au(W0) [121 H51 H3| H5] . 

Square lattices with periodical boundary conditions are used in the present 
work to represent the surface of the electrode. Each lattice point represents 
an adsorption site for an adsorbate or a substrate atom. The former may 
adsorb or desorb on each randomly selected site, while the latter may only 
move on the surface jumping from the selected site to a neighboring one. In 
this way, our model corresponds to an open system for one of its components, 
that is, the adsorbate. This has physical correspondence with the setup of 
the electrochemical experiment, where only the metal in equilibrium with 
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its cations in solution may dissolve or be deposited in the potential range 
considered. 

Concerning the adatoms of the same kind of the substrate, some consider- 
ations must be made regarding the existence of surface defects. These atoms 
may in principle move freely on the surface and minimize the free energy 
of the system by a number of mechanisms. For example, isolated substrate 
atoms may heal defects through their incorporation to a defective cluster, or 
small substrate islands may dissolve to join large ones, like shown in previous 
simulation work by Stimming and Schmickler |16| . Thus, different substrate 
structures can be imposed as initial conditions for each simulation. In the 
present case, we employ islands of substrate atoms on the surface obtained 
by means of simulated annealing techniques [9]. This was undertaken in or- 
der to emulate some of the defects that can be found on a real single crystal 
surfaces, like kink sites, steps, etc. 

On purely energetic grounds, metal electrodeposition on clean metal sur- 
faces with islands should take place according to the following sequence: 

i) Decoration of the lower part of island edges. This is so because the 
binding energy is the lower (more negative), the higher the coordination of 
the site where a metal atom becomes adsorbed. 

ii) Deposition on the rest of the flat surface around the islands. 

iii) Some nucleation should occur at the top of the islands. This requires 
a more negative potential than that for island decoration (i) or flat surface 
filling(ii) because the border of the growing nuclei makes them less stable 
than a growing monolayer (nucleation overpotential) 

iv) The border of the islands at the top should be finally decorated. These 
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sites are particularly unstable due to their low coordination. 

In the present work we are mainly interested in island decoration and 
the possible mixing between island atoms and depositing atoms following 
immediately after this decoration. This involves steps i) and ii) described 
above. Steps iii) and iv) could be of course simulated but this will not be 
considered in the present work. 

2.2 Energy calculations 

The energy E of the system is related with the classical lattice-gas model 
with pair potential interactions between nearest neighbors: 

+ J226 

(a,/?) 

(1) 

where the symbol S represents the Kronecker delta and c a indicates the oc- 
cupation number of the a site [corresponding to a pair of (x,y) coordinates]. 
Unlike the classical lattice-gas model, the occupation number may assume 
three different values: 0, 1 or 2, if the corresponding site is empty, occupied 
by an atom of the same kind of the substrate or occupied by an adsorbate 
atom, respectively. The notation (a, (3) implies a sum taking account of all 
pairs of nearest neighbors, is the adsorption energy of a particle of type % 
and Jij are the pairwise interactions between nearest neighbors. 

The values of the parameters considered here for the different systems are 
summarized in table 1. These values were obtained according to calculations 
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employing the embedded atom model [El El HO] . In a previous work [10], we 
have calculated the adsorption energies corresponding to the different envi- 
ronments of the adsorption site, taking into account first, second and third 
neighbors. Some of the possible configurations were ilustrated in Fig. 1 and 
corresponding energies were indicated in table 1 of reference pU]. The values 
employed in the present work for Cj were based on the adsorption energies 
corresponding to configuration 1 (that is, adsorption on a site without neigh- 
bors). The values employed for the parameters were based on the others 
configurations (2, 3, etc., which correspond to different occupation of neigh- 
boring sites) assuming pair potential interactions between nearest neighbors 
and taking an average value. 

2.3 Grand Canonical Monte Carlo 

Square lattices of M = L x L adsorption sites with periodical boundary 
conditions are used here to represent the surface (in this case, L = 100). In 
the case of adsorption isotherms on clean surfaces, the initial state consists of 
2000 adsorbate atoms distributed at random (that makes a coverage degree 
of 0.2). For the case of adsorption in the presence of surface defects, the 
initial state is characterized by Nj, = 1000 atoms of the same nature as that 
of the substrate forming islands previously generated by simulated annealing 
[HI HO], and 2000 adsorbate atoms distributed at random. The islands of 
substrate atoms are useful to emulate surface defects like kink and step sites. 

The simulation consists in the realization of a certain number of Monte 
Carlo Steps (MCS) in order to equilibrate the system and then another set 
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of MCS evaluation, from time to time, the thermodynamic quantities of 
interest (like energy of the system or coverage degree) in order to obtain the 
average value. This is performed for fixed values of temperature and chemical 
potential (Grand Canonical ensemble). 

Each MCS implies the realization of M trials. Each trial consists on the 
random selection of one lattice site, with the realization of one of the three 
following processes: 

a) If the site is empty (occupation 0), the creation of an adsorbate atom 
is attempted (yielding occupation 2). 

b) If the site is occupied by an adsorbate atom (occupation 2), its desorp- 
tion is attempted (occupation 0). 

c) If the site is occupied by a substrate atom (occupation 1) a nearest 
neighbor site is selected at random. If this new site is empty, a move 
to this position is tried. 

In the three cases, the change is accepted with probability: 



AE - fiAN a 

P = mm{l, exp[ — ]} (2) 



where AE (AN a ) represents the difference between the energies (number 
of adsorbed particles) of the final and initial states, ks is the Boltzmann 
constant, T is the temperature and \i is the chemical potential. Note that 
the value of AN a is +1 in case a), —1 in case b) and in case c). 
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The algorithm described above has some differences with the algorithm 
employed in Refs. (9j [JO]. In the present case, only the substrate atoms 
can move between adjacent sites, although this possibility is indirectly con- 
templated for adsorbate atoms in the case of desorption from one site and 
adsorption in a neighboring site. This difference is not essential, the main 
target was to save computational time. 

2.4 Some definitions 
2.4.1 Coverage degree 

The coverage degree, #(//), which is defined as the number of sites occupied by 
an adsorbate atom divided by the number of available sites (all adsorption 
sites that are not occupied by a substrate atom), is obtained as a simple 
average: 

where (N a ) is the mean number of adsorbate atoms on the surface and (• ■ ■) 
denotes the time average over the Monte Carlo simulation runs. 

In the case of adsorption in the presence of surface defects, that is, when 
some substrate islands are present, we also define some partial coverages, like 
the coverage degree of steps sites, 8 S , and the coverage degree of kink sites, 
9k- 



<?U (4) 
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where (iV S)t ) ( (N k;t ) ) is the average number of step (kink) sites on the surface 
and (N Sj0 ) ( (N^ ) ) is the average number of step (kink) sites occupied by 
an adsorbate atom. 

2.4.2 Differential heat 

We calculate the differential heat of adsorption, defined as [T8] : 



9(E) 

qd = ~dm (6) 

where (E) is the mean energy of the adlayer. 

The physical interpretation of this quantity is the energy associated with 
removing a particle from the surface at each moment. 

In the case of computer simulations, the most appropriate way of per- 
forming the calculation of this value is by means of averaged quantities. We 
employ the following formulation [TT1 19|: 

(EN a )-(E)(N a ) 
2.4.3 Quantification of island disintegration 

As discussed below, for some systems the substrate islands are found to 
disintegrate upon adatom deposition. In order to quantify this phenomenon, 
we define some (normalized) quantities characterizing the disintegration of 
substrate islands present on the surface when the adsorbate particles adsorb, 
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as follows: 



• Fraction number of 1 — pairs: 



MP E<g,/?>fl(l,c/?) x (5(0, c a ) 



where the sum yields the number of adjacent substrate atom-empty 
site pairs, c is the lattice connectivity (in this case c = 4). 

Fraction number of 1 — 2 pairs: 

JVPl - 2 " dv^ (9) 

where the sum yields the number of adjacent substrate-adsorbate pairs. 
Total number of 1 — X pairs, with X = 0, 2, that we define from: 



NP T = NP ± _ + NP 1 _ 2 (10) 

which corresponds to the fraction of sites surrounding substrate atoms 
that are not occupied by particles of the same kind (type 1). The re- 
sulting quantity, further on called "fraction of island disintegration", is 
given by the number of substrate atoms/adsorbate atoms and substrate 
atoms/empty sites pairs. 

The maximum value of NPt will be reached in the case where the islands 
are completely disintegrated, that is, when there is a minimum of neighboring 
substrate-substrate pairs. In principle, the upper value of NPt could be 
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1, in the case where no substrate-substrate pairs occur at the maximum 
coverage by the adsorbate. However, if we only take into account statistical 
considerations, it is to be expected that the probability of finding a neighbor 
of the same species (1) should be given by the fraction N d /M 7 which in 
the present simulation is fixed to the value equal to 0.1. In other words, 
it is expected that in the limit of high temperatures where entropic effects 
prevail, and for large 9, we should get that NPt — ► 0.9. 

For a given distribution of well equilibrated substrate islands, the mini- 
mum value of NPt will be given by the fraction of pairs involving atoms that 
belong to the edge of the island at beginning of the simulation. If the islands 
do not disintegrate at all, remaining intact, NP T should remain the same 
all over the simulation. On the other hand, if the islands disintegrate, NPt 
should reflect the increasing connectivity of the substrate atoms to sites of a 
diferent kind. 

3 Results and discussion 
3.1 Adsorption on clean surfaces 

Adsorption isotherms were simulated for the systems Ag/Au(100), Ag / Pt(100) , 
Auj Ag (100) and Pt/Ag (100). The corresponding curves are plotted in the 
upper part of Figs. 1 and 2 (a) and b)), for different temperatures. It can be 
seen that at low temperatures the isotherms show an abrupt jump, typical 
of first order phase transitions. As the temperature increases, the isotherms 
become smoother, specially for the first two systems. 
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The present model is analogous to the Ising model in the sense that lateral 
interaction between adsorbing particles is considered as a pair-potential and 
only between nearest neighbors. It is well known that in this case the critical 
temperature for 9 = 0.5 can be estimated as |20j: 



2k B ln{V2-l) 

where J is the lateral interaction between adsorbing particles. 

Taking into account the values of J 22 employed here (see table 1), the 
estimated critical temperatures are: Tc = 1843K for the system Ag/Au(100); 
T c = 1382iT for Ag/Pt(100); T c = 3028K for Au/Ag (100) and T c = 3950X 
for Pt/Ag (100). This is in agreement with the appreciation of the isotherms 
in Figs. 1 and 2, where it can be seen that the critical temperature must 
be between lOOO-RT and 2000^ for the systems Ag/Au(100) and Ag/Pt(100), 
close to 3000fr for Au/Ag {100) and above 300CLK" for Pt/Ag {100). 

The chemical potential at which the isotherms of figure 1 and 2 present 
the step, say /is, provides a measure for the affinity of the adsorbate for the 
substrate. This quantity can be substracted from the corresponding binding 
energy of the adsorbate, and compared with the so-called underpotential shift 
A(/) up d, which was first defined by Kolb et al. [1] as the potential difference 
between the desorption peak of a monolayer of a metal M adsorbed on a for- 
eign substrate S and the current peak corresponding to the dissolution of the 
bulk metal M. The magnitude of A0 up d is a measure of the affinity of the 
adsorbate for the substrate, as compared with the affinity of the adsorbate 
with itself, and can be written in terms of the chemical potential per parti- 
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cle of the atom adsorbed on a foreign substrate, /x[(5)M],and the chemical 
potential of the same species in the bulk fi[(M)M], according to|21J: 

A0 upd = — (/i[(M)M] - fi[(S)M]) (12) 

In the present model, fi[(M)M] can be replaced by the cohesive energy of 
the adsorbate, and fi[(S)M] may be in turn substituted by the chemical 
potential at which the isotherm presents the step /x^, since this chemical 
potential correspond to the occurrence of the adsorbate phase. 

A0u P d(^l) = -^(^r-^]) (13) 

2:60 

where we have introduced the superscript T to denote that this is a the- 
oretical prediction for the (lxl) adsorbate structure . On the other hand, 
the experimental estimation of the underpotential shift A0^ d can be made 
from: 



ACpd(l^) = 4(1X1) - 4>Nernst (14) 

where 0(1x1) denotes the potential at which the (lxl) adsorbate phase ap- 
pears, and 4>Nemst is the reversible deposition potential for the electrolyte 
solution employed. The Ag/Au(100) system has been considered by dif- 
ferent groups [121 CEU E2]. Ikemiya et al reported (pNemst = —58 mV and 
0(1x1) ~ 180m V vs Hg / Hg 2 SO l i for this sytem, thus yielding A0^ d (lxl) ~ 
238 mV . In the case of the Ag/Pt(100)system, Aberdam et al.[23] reported 
that the peak for the deposition of the first monolayer of Ag was located 0.48 
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V more positive than the Nernst reversible potential, so that we take this 
value for A.<p^ d (lxl). We show in Table II these experimental values along 
with the present theoretical estimations and those of our previous work using 
the more accurate embedded atom method potentials. It can be appreciated 
that the pair potential approximation performs almost as well as the more 
sophisticated many body one. There is no evidence in the literature for un- 
derpotential deposition of Au on Ag(100) or for Pt underpotential deposition 
on Ag(100). A related system, Pt deposition on Au(100) has been measured 
by Waibel et al.[24], with the finding that Pt deposition takes place at over- 
potentials. Since the cohesive energy of Au is larger than that of Ag, it is 
expected that if Pt upd is not found on Au(100), it will be even less probable 
on Ag(100). 

The lower parts of Figs. 1 and 2 (parts c) and d)) show the differential 
adsorption heats. The qualitative form of the curves of qd vs. /i is similar to 
the adsorption isotherms (9 vs. /i). The observed values of qd are in agree- 
ment with the estimation of qd ~ —62 at very low coverages (corresponding to 
the adsorption of particles in sites without neighbors) and qd ~ —£2 — 4 x J 2 2 
at high coverages (corresponding to the adsorption of particles in sites sur- 
rounded by four nearest neighbors of the same kind). 

3.2 Adsorption in the presence of substrate islands 

Figs. 3 and 5 show the adsorption isotherms for two of the four studied 
systems in the presence of surface defects (substrate atoms present in the 
monolayer forming islands). In the figures, the isotherms for the complete 
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monolayer are plotted together with the isotherms of steps and kink sites for 
three different temperatures [parts a), b) and c)], as well as the differential 
adsorption heats for the three temperatures [part of)]. 

It can be seen that for Ag/Pt(100) (Fig. 3) and Ag/Au(100) (not shown), 
the kink sites are occupied first, then the step sites and finally the complete 
monolayer. For the other two systems, at T = 300-fT, there is no differentia- 
tion between the different kind of sites, but at greater temperatures a small 
difference arises (see parts a) — c) of Fig. 5 for the system Au/Ag(100), the 
system Pt/Ag(100) is not shown, but the results are very similar). 

Analyzing the behavior of the differential heat of adsorption [part of)], it 
can be seen that for the first two systems (see Fig. 3 for Ag/Pt(100)), at 
low temperatures, there are four differentiated stages, before the adsorption 
of the complete monolayer (last part). The first stage with qd ps — (e 2 + 2J i2 ) 
(filling of substrate kink sites); the second one with qd ~ — (e 2 + Ji 2 + J 22 ) 
(filling of step sites besides to another adsorbate particle); the third one with 
q d ps — (e 2 + 2J 22 ) (filling of an adsorbate kink site, that is, adsorption on 
terrace sites next to two adsorbate particles during monolayer growth) and 
the fourth one with qd ~ — (e 2 + 4J 22 ) (monolayer completion). All these 
stages are illustrated in Fig. 4. These estimations are closest to the results 
of the simulations for the system Ag/Pt(100). 

For the other two systems (see Fig. 5 for Auj Ag (100)) there are at first 
sight only two more or less well defined stages: before and after monolayer 
completion. In the first part, where only a few adsorbate particles come into 
the system, the differential heat can be written as qd ~ — [e 2 + 4(J 12 — J n ) + 
2Jn]. This value corresponds to the replacement of a substrate particle by 
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an adsorbate particle in the middle of an island and the positioning of the 
former at the edge of the island. The interpretation of the value of qa is not 
straightforward for monolayer completion. 

Fig. 6 shows the final state of a portion of the simulation cell at three 
different chemical potentials for the system Ag/Au(100) in the presence of 
defects at two different temperatures. At T = 300i^, the islands remain 
almost unchanged upon adsorbate deposition and the sites are filled follow- 
ing the order: 1) kinks; 2) steps; 3) terraces. At T = 1000-fT, the general 
tendency remains but the islands show a certain disintegration. The same 
behavior is found for the system Ag / Pt(100) (not shown). 

On the other hand, in the case of the system Au/ Ag(100) at T = 300K 
(upper part of Fig. 7), the islands do not change their shape significantly 
but some adsorbate atoms penetrate inside them. At T = 1000-fT (lower part 
of Fig. 7) the islands disintegrate completely. The same general behavior is 
found for the system Pt/Ag (100) (not shown). 

Figs. 8 and 9 show the adsorption isotherms and the fraction of island 
disintegration (characterized by the value NP Tl defined in equation [TOl) in 
the presence of substrate islands at three different temperatures for two of 
the considered systems. 

In the case of the system Ag/Au(100) (Fig. 8), it can be verified that, 
at low temperatures, the islands remain unchanged, since the value of NP? 
is constant along the whole isotherm. As expected, the value of NPi_ 
decreases and the value of iVP]_ 2 increases as the adsorbate atoms cover the 
edges of the islands. At higher temperatures, the value of NP T increases when 
the adsorbate atoms enter. That means that both atom classes mix with 
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each other, that is, the substrate islands disintegrate. At an intermediate 
temperature (T = lOOOif) the value of NPt presents a maximum close to 
the critical chemical potential due to the greater mobility of the atoms in 
that situation. 

For the system Ag / Pt(100) (not shown) the general trend is very similar 
to the previous case. 

The adsorption isotherms for the systems Au/Ag (100) (not shown) and 
Pt/Ag (100) at different temperatures (Fig. 9a) show that, as the tempera- 
ture increases, the jump in coverage moves towards more negative chemical 
potentials. Considering the fraction of island disintegration, an important 
difference with respect to the other systems consists in that even at low tem- 
peratures the islands start to disintegrate when the adatoms are deposited 
(Fig. 96). This is evident in the increase of the quantity NPt- Another 
important detail to be emphasized is that for high temperatures island dis- 
integration is considerable, even for low adsorbate coverages (Fig. 9d). As 
anticipated in section 2.4.3, it can be observed in Figs. 8d and 9d that for 
all systems NP T — > 0.9 in the limit of high temperatures and adsorbate 
coverage degrees. 

The general picture that we get from the present results is that as long 
as the susbtrate has a binding energy that is considerably larger than that 
of the adsorbate, the islands of the former remain relatively unaltered upon 
adsorbate formation. As the binding energy of the substrate approaches 
that of the adsorbate or becomes smaller, the islands become unstable and 
disintegrate relatively easily. The study of Pd (E coh = —3.91 eV) deposition 
on Au (100) (E^h = —3.93 eV)shows in fact the formation of a surface alloy in 
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the underpotential region [25]. The study of Pt (E coh = —5.77 eV) deposition 
on Au(100) shows some remarkable features. Although surface alloying was 
not reported for this system [21], the Pt deposit exhibit irregular shapes and 
recent molecular dynamic simulations point towards the existence of a surface 
alloy [26]. 

4 Conclusions 

In the present paper we revisit sytems previously simulated in Refs. P, [TO] 
with a different model for the metal-metal interactions. Instead of the com- 
putationally demanding many-body potentials employed there, we use here 
attractive pairwise aditive potential interactions between nearest neighbors 
with only a minimal set of parameters. While the main results remain quali- 
tatively similar, the present modelling is considerably simple and workable in 
the framework of the pairwise additive potential lattice model. It also saves 
considerable computational time. Furthermore, we have considered here the 
influence of the temperature and some new quantities, like differential heat 
and fraction of island disintegration, that helped us in the understanding of 
the simulations. 

The present formulation shows that the main characteristics to be taken 
into account is the strenght of the interactions between adsorbate atoms as 
compared with the interaction between substrate atoms. In this sense, we 
can divide the four systems studied as examples into two groups. The first 
group is integrated by the systems Ag/Au(100) and Ag/Pt(100). In these 
cases the interactions between adsorbate atoms is weaker than the interaction 
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between substrate atoms (i.e.,| J22I < |^n|)- The second group is integrated 
by the systems Au/Ag (100) and Pt/Ag (100), where the interactions between 
adsorbate atoms is stronger than the interaction between substrate atoms 
(i.e., |J 22 | > |Jn|). 

For defect-free surfaces, the adsorption isotherms show an abrupt jump 
at low temperatures and become smoother at high temperatures. This is 
indicative of the existence of a first order phase transition and a critical 
temperature related with the interaction between adsorbate atoms. 

The adsorption in the presence of islands of the same nature as that of 
the substrate on the surface was also studied. 

For systems in which the interaction between adsorbate particles is weaker 
than the interaction between substrate particles (Ag / Au(100) and Ag / Pt(100)) , 
the substrate islands remain relatively unchanged at low temperatures and 
show a certain degree of mobility at high temperatures. For systems in which 
the interaction between adsorbate particles is stronger than the interaction 
between substrate particles (Au/Ag (100) and Pt/Ag (100)), the adsorbate 
atoms penetrate into the islands at low temperatures and the islands are 
completely disintegrated at high temperatures. 

The simplicity of the present formulation will allow the analysis of prob- 
lem at hand in terms of a few parameters than can be systematically varied. 
In other words, instead of system-oriented simulations, studies with these 
control parameters can be performed, with the consequent gain of generality. 
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6 Figure Captions 

Figure 01 a) Adsorption isotherms at different temperatures for the system 
Ag/Au(100). b) Adsorption isotherms at different temperatures for the sys- 
tem Ag/Pt(100). c) Differential heats for the adsorption isotherms of the 
system Ag/Au(100). d) Differential heats for the adsorption isotherms of the 
system Ag/Pt(100). 

Figure 02 a) Adsorption isotherms at different temperatures for the 
system Au/Ag (100). b) Adsorption isotherms at different temperatures for 
the system Pt/Ag (100). c) Differential heats for the adsorption isotherms of 
the system Au/Ag (100). d) Differential heats for the adsorption isotherms 
of the system Pt/Ag(100). 

Figure 03 a) — c) Adsorption isotherms in the presence of surface defects 
for the complete monolayer, the step sites and the kink sites at three different 
temperatures for the system Ag/Pt(100). d): differential heat of adsorption 
at the three temperatures for the same system. 

Figure 04 Schematic representation of the top view of four environment 
types close to the adatom adsorption site, a): substrate kink site, b): sub- 
strate step site, next to another adsorbate particle, c): adsorbate kink site. 
d): hollow site (upon monolayer completion). Filled circles denote substrate 
atoms. Empty circles represent adsorbate atoms. The empty dashed circle 
denotes the adsorption site under consideration. The underlying substrate 
atoms are not shown. 

Figure 05 a) — c)Adsorption isotherms in the presence of surface defects 
for the complete monolayer, the step sites and the kink sites at three different 
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temperatures for the system Au/Ag (100). d): differential heat of adsorption 
at the three temperatures for the same system. 

Figure 06 Snapshots showing the final state of the surface at different 
chemical potentials for the system Ag/Au(100) in the presence of surface 
defects at T = 300-K" and T = lOOO-ft'. Filled circles represent gold atoms 
while non-filled ones represent silver atoms. 

Figure 07 Snapshots showing the final state of the surface at different 
chemical potentials for the system Au/Ag (100) in the presence of surface 
defects at T = 300K and T = lOOOi^. Filled circles represent silver atoms 
while non-filled ones represent gold atoms. 

Figure 08 . a) Adsorption isotherms in the presence of surface defects 
for the complete monolayer at three different temperatures for the system 
Ag/Au(100). b) - d) Fraction number of pairs (NP^o, NP^ 2 and NP T ) at 
the same temperatures. 

Figure 09 a): Adsorption isotherms in the presence of surface defects 
for the complete monolayer at three different temperatures for the system 
Pt/Ag(100). b) - d) Fraction number of pairs (NP^ , NP X - 2 and NP T ) at 
the same temperatures. 
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7 Tables 



System 




£2 


Jn 


J22 


J12 


Ag/Au(lOO) 


-3.05 


-2.58 


-0.54 


-0.28 


-0.42 


Ag/Pt(lOO) 


-4.34 


-3.13 


-0.83 


-0.21 


-0.56 


Au/Ag (100) 


-2.39 


-3.11 


-0.25 


-0.46 


-0.36 


Pt/Ag(lOO) 


-2.30 


-4.22 


-0.17 


-0.60 


-0.41 



Table 1: Parameters representing the adsorption energies and the inter- 
action energies (in eV units) employed here for the considered systems. 



System A«^ pd (lxl) / V A«/£ pd (lxl), EAM / V A^ A (lxl) / V 

Ag/Au(100) 0.23 0.17 0.24 a 

Ag/Pt{100) 0.67 0.55 0.48 b 

Au/Ag(100) 0.04 -0.08 < 0? c 

Pt/Ag(100) -0.44 -0.53 < 0? c 

Table 2: Calculated and experimental underpotential shifts for the sys- 
tems considered in the present work. A0^ pd (lxl) denote results of this 
work, A0^ pd (la;l) / EAM are results from previous simulations using the 
many-body potential from the embedded atom methodfTT]. A0^ p p d (la;l) are 
experimental estimations taken from the literature. A0^ d (lxl) values were 
calculated using equation ([H]), with E%> h = -2.85, -3.93 and -5.77 eV for Ag, 
Au and Pt respectively. 

a) taken from reference [22], b) Taken from reference [23] c)no upd has 
been reported in the literature for these systems so far. 
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